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Abstract. We present evidence for the existence of shock 
waves caused by the formation of the large-scale structure. 
In some clusters of galaxies peripherally located sources of 
extended diffuse radio emission exist, the so-called cluster 
radio relics. They have steep radio spectra but no appar- 
ent cutoff, as old remnants of radio galaxies usually have. 
Therefore particle acceleration has to take place within 
them. We propose that shock structures of the cosmolog- 
ical large-scale matter flows are responsible for the accel- 
eration of relativistic electrons: cluster accretion shocks 
and bow shocks of merger events. We develop a theory 
of radio plasma having traversed these shocks and com- 
pare it to observational data of nine radio relics (0038-096, 
0917-h75, 1140-1-203, 1253-1-275, 1712-^64, 1706-1-78, 2006- 
56, 2010-57, 1401-33) and their host clusters (A85, A786, 
A1367, Coma, A2255, A2256, A3667, S753). The neces- 
sary accretion power, the spectral index of the radio spec- 
trum, the acceleration efficiency of the shock, the diffusion 
coefficient in the post-shock region, and the predicted ra- 
dio polarization in all of our examples fit into a coherent 
interpretation of the observational data. Since polariza- 
tion measurements are available only for four sources, the 
predictions of our theory can be independently checked 
using other examples. The predicted values of the shock 
compression ratio, density and temperature of the infalling 
gas, magnetic field strength of the shocked and unshocked 
radio plasma are discussed within the frame of structure 
formation theory. 
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1. Introduction 

Large extended sources of synchrotron emission are found 
in some clusters of galaxies: radio halos and radio relics 
(Fig. ^ for reviews, see Jaffe 1992; Feretti & Giovannini 
1996). The rare cluster radio halo phenomena are large 
regions of diffuse unpolarized radio emission, centrally lo- 
cated, and have radio structures roughly similar to that 
of the thermal X-ray emission (Deiss et al. 1997). The so- 
called radio relics are peripherally located, and are more 
irregularly shaped polarized radio sources. Radio relics are 
believed to be the remnants of radio lobes of radio galax- 
ies, where the former active galaxy has become inactive 
or has moved away. 

For radio halos a number of candidate sources for the 
relativistic, synchrotron emitting electrons were discussed 
in the literature (Schlickeiser et al. 1987 and references 
therein). Yet, the source of the energetic electrons in the 
relics is unclear, due to the shortness of the electron cool- 
ing time compared to the time since electron injection 
from a galaxy at the location of the relic was possible. 
Since the radio spectra of relics are steep but frequently 
do not exhibit any cutoff in the observed range, an effi- 
cient particle acceleration mechanism has to be present 
within or close to them. We show that shocks of the large- 
scale gas motion are expected at the typical peripheral 
locations of relics, either resulting from cluster mergers, 
or from steady state accretion shocks by gas falling into 
the cluster potential. Kang et al. (1997) showed that pro- 
tons might be accelerated in cluster accretion shocks to 
energies comparable with the most energetic cosmic ray 
events observed at earth. Also electrons should be accel- 
erated and become visible by radio emission at locations 
where magnetic fields are present. Therefore the magne- 
tized plasma accumulated behind the shock or left behind 
as a remnant of a radio galaxy, could furnish this accel- 
eration region and should also increase the efficiency of 
the acceleration mechanism. We demonstrate that for the 
case of the radio relic 1253-1-275 close to Coma A (Fig. |]), 
the observed spectral index and polarization of the radio 
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Fig. 1. Westerbork Synthesis Radio Telescope map at 327 MHz of the Coma cluster from Giovannini et al. (1991). The central 
halo source Coma C and the cluster relic 1253+275 are identified. Capital letters indicate some extended Coma cluster galaxies, 
as following: A = NGC 4839, B = NGC 4827, C = NGC 4789. 10' corresponds to 400 kpc /I'V 

These estimates depend on the assumption of a spher- 
ical flow, of a polytropic index of the gas of 7gas = 5/3, 
and they depend weakly on the assumed Einstein-de Sitter 
cosmology. Since the real infall pattern to clusters should 
be aspherical and correlated with the surrounding large- 
scale structure (Colberg et al. 1997), the shock radii and 
the shock velocities are expected to differ from these es- 
timates. They might be regarded as average quantities, 
where from a single cluster can easily deviate. Inserting 
the X-ray temperature of the Coma cluster of 8.2 keV 
(Briel et al. 1992) and the redshift of z = 0.0233 (Fadda 
et al. 1996) gives a shock radius of = i.Sh^^Mpc and 
a shock velocity of 2.0 ■ lO^kms^^. If 1253-1-275 is at the 
position of this shock, the projected distance to the cluster 
center of 2.9/i5o^Mpc (Giovannini et al. 1991) implies an 
angle between line-of-sight and normal of the shock front 
of 5 = 37°. The true angle is ambiguous (180° — 6 or S), 
and depends on the location of the relic on the line-of-sight 
(foreground or background with respect to the cluster cen- 
ter). A possible configuration is sketched in Fig. ^ , where 
the relic is located on the back side of the cluster accre- 
tion shock, since then a physical connection to the nearby 
radio galaxy NGC 4789 is possible as explained in Sect. |[ 



emission agree well with the prediction for an accretion 
shock. Eight other relics are also briefly discussed. 

We adopt for the following fl^ = 1 and Ho — 
50 /i5o kms~iMpc~^ 

2. Particle acceleration at cluster accretion shocks 
2.1. Geometry 

Properties of the accretion shocks of clusters of galaxies, 
the shock radius rg and the velocity Vg of the infalling mat- 
ter measured in the shock frame, are related to the depth 
of the gravitational potential of the cluster. Since the clus- 
ter gas temperature also depends on the potential, the 
shock properties can be expressed in terms of the observed 
temperature fcTobs- Using results of one-dimensional sim- 
ulations of accreting flows onto clusters (Ryu & Kang 
1997), Kang et al. (1997) derived 
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Fig. 2. Sketch of the assumed geometry in the Coma cluster. 
The accretion flow is indicated by arrows. The shaded region is 
the intra-cluster medium. The dotted area represents regions 
of radio plasma. The line-of-sight and the line connecting the 
relic and the cluster center are lying within the plane of the 
figure. NGC4789 is placed upstream above the radio relic as 
a source of relativistic electrons, as discussed in Sect. ^. Its 
line-of-sight velocity component points away from the observer. 
All positions could be different due to shifts parallel to the 
line-of-sight. 

The length of the relic is 1.2Mpcft,5Q^, whereas the thick- 
ness cannot be given directly from radio observations. We 
assume a thickness of I? = 0.1 Mpc/i^Q^, which is smaller 
than the observed projected thickness due to the action of 
shock compression in the infall direction. A detailed model 
of the structure of the relic would allow to deproject the 
apparent surface. For only moderately large projection an- 
gles {6 < 60°) resulting deprojected surfaces are not much 
bigger than the observed ones. Since the geometric quan- 
tities are only used for order of magnitude calculations of 
the infalling kinetic energy flux, the shock efficiency, and 
the diffusion coefficient within the radio plasma, we pre- 
fer not to specify a detailed geometrical model. We equate 
projected to real surface of the relic, which is a weak un- 
derestimate. The relic is slightly V-shaped and the radio 
emission is more sharply edged on the outer side (Gio- 
vannini et al. 1991). The V-shape could result from an 
aspherical shock, but would also be easily explained by an 
intrinsic shape of the magnetized plasma, since the relic 
is not seen edge-on. The edges of the observed projected 
radio structure should differ. The radio plasma seen at the 



inner edge (left edge in Fig. ^ is more distant to the shock 
than that at the outer edge (right). The outer edge should 
be sharper due to the better confined plasma at the shock 
side, and it should have a lower than average spectral in- 
dex, since the electrons there are reaccelerated recently. 
Whereas the inner edge is expected to be smoother and 
to have a higher than average spectral index. The reason 
for the latter is that the reaccelerated electrons within the 
radio plasma seen there, at a position more distant to the 
shock, had more time for cooling. 

2.2. Compression ratio 

The radio spectral index a {Si, ^ of 1253-1-275 is 

a = 1.18 ± 0.06 (Giovannini et al. 1991) and is therefore 
larger than a = 1, expected for a strong shock accelerat- 
ing electrons, with simultaneous synchrotron and inverse- 
Compton losses. Larger spectral indices arc found at the 
inner edge of the relic, showing evidence for steepening 
of the electron momentum distribution due to the higher 
age of the electrons seen there, more distant to the accel- 
eration region of the shock. The radius of the accretion 
shock structure (assumed to have a spherical surface) is 
large compared to every intrinsic length scale of the accel- 
eration process. Thus the theory of planar shocks can be 
applied. The spectral index can be explained by allowing 
the shock compression ratio R to be lower than the canon- 
ical value i? = 4 of a strong shock of gas with a polytropic 
index of 7gas = 5/3. The radio spectral index a of an equi- 
librium electron population accelerated and cooled at the 
same time, is related to the shock compression ratio by 
(e.g. Drury 1983): 



where a polytropic index of 7gas = 5/3 is assumed. The 
compression ratio is therefore R = 3.2 ± 0.2. From the 
theory of shocks (Landau & Lifschitz 1966) the pressure 
and temperature ratio between the down- and up-stream 
region (inside and outside the cluster shock front, in the 
following denoted by P2 and Pi, etc.) can be derived: 

^ 4i?-l ^ a+l Tl^LEl U\ 
Pi 4-R a - 1 ' Ti R Pi ' ^ ' 

For 1253+275 the resulting pressure ratio is P2/P1 = 
11 - 22. The temperature ratio T2/T1 = 3.8 - 6.4 re- 
quires that the infalling matter has a temperature of 
Ti ^ 0.6 — 1.1 keV, assuming the post-shock temperature 
to be roughly half the central cluster temperature. Simula- 
tions predict this radial temperature decrease (Navarro et 
al. 1995), and observations indicate a radially falling tem- 
perature profile in clusters (Honda et al. 1996; Markevitch 
1996; Markevitch et al. 1996, 1998). 

This temperature of the infalling matter is reason- 
able, since this gas should flow mainly out of sheets and 
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filaments of the cosmological large-scale structure, and 

therefore was preheated by the accretion shocks at the 
boundaries of these structures. A simulation of cosmolog- 
ical structure formation by Kang et al. (1996) shows that 
typical temperatures of filaments are above 0.1 keV. Adi- 
abatic compression and internal shocks within the flows 
onto clusters might raise these temperatures to a level of 
0.5 — 1 keV, which is needed in order to explain the steep- 
ness of the synchrotron spectrum of cluster relic sources 
caused by weak shocks. 

The temperature of the infalling matter gives the 
sound velocity ci, which enters into a second estimate of 
the shock velocity, using the theory of planar shock waves 
(Landau & Lifschitz 1966): 

I (7gas - 1) + (7gas + l)P2/Pl ... 
I's, predicted — ^1 \ • [0) 

This estimate gives a shock velocity for Coma of 
Ks. predicted = 1850 km/s, slightly smaller than that fol- 
lowing from the theory of Kang et al. (1997). The differ- 
ence vanishes if a smaller temperature drop of w 1.6 from 
the center to the shock radius is used instead of a factor 
of 2. Comparing both estimates shows self-consistency of 
our model, but this is not a completely independent test, 
since both values depend mainly on the same quantity 

(y, ~ (fcTobs)^/^). 

In the model of acceleration of ultra-high-energy cos- 
mic rays at cluster accretion shocks of Kang et al. (1997) 
a compression ratio close to i? = 4 is favored in one partic- 
ular model to fit the shape of the spectrum of ultra-high- 
energy cosmic rays using a low and energy independent 
escape efficiency of cosmic rays from the cluster against 
the upstream flow. However, a lower compression ratio 
(and therefore a steeper cosmic ray production spectrum) 
can be compensated by taking a reasonable energy de- 
pendence of the escape efficiency into account. In fact, in 
Kang et al. (1997) the acceleration and escape efficiency 
was rather low, with w lO^"' giving a reasonable fit only to 
the observed high-energy cosmic ray data. Using the com- 
pression ratio of i? = 3.2 ± 0.2 gives a proton spectrum of 
^-2.36±o.i2 ijnpiying for the same flux at a few 10^^ eV a 
combined efficiency which is in the range 0.03 — 0.1 (Fig. 
5 in Kang et al. 1997) and so consistent with estimates for 
supernova remnants (e.g. Drury et al. 1989; Volk 1997). 

2.3. Shock Efficiency 

A rough estimate of the gas densities on both sides of 
the shock sphere can be gained by extrapolating the 
electron density profile of the Coma cluster, viz. rie = 
no,o [1 + {r/rcoT:c)^]~^^^'^, with the parameters ne,o = 
3 • 10-3/i5Q^/^cm-3, rcore = 400 /j-^ kpc, and P = 0.75 
(Briel et al. 1992). This gives a gas density of the order 
ne,2 ~ 10~^ h^Q^^cm.~^ at a cluster radius of 4.8 h^Q Mpc. 
The density of the infalling matter should therefore be of 



the order of ne,i « 3.5-10"^ h^Q^^cm.~^. This is consistent 
with the somewhat lower gas density in the intergalactic 
space of < 10~^ cm~^ far away from clusters, derived e.g. 
from radio observation of the giant radio galaxy 1358-1-305 
(Parma et al. 1996). 

The projected extent of 1253-^275 is 1.2 x 0.4 h^^ Mpc^ 
(Giovannini et al. 1991), but its real surface S exposed to 
the accretion flow might be larger. The kinetic power of 
the infalling matter on the relic 

Qflow = \ ne,i mp y/Fs 5 « 5 • 10^^ h^^^^ erg s'^ , (6) 

is mainly transformed into thermal energy, but some frac- 
tion is convert(xl into relativistic particles. Vg — Vg {R — 
I) / R is the velocity of the infalling matter measured in 
the cluster's inertial frame. The integrated radio power 
between 10 MHz and 10 GHz of 8 • 10""' h^^ ergs^^ (Gio- 
vannini et al. 1991) is three orders of magnitude lower, 
and can be easily powered by the dissipative processes in 
the shock. Assuming a field strength of 1/uG implies that 
only 10% of the electron radiation losses are synchrotron 
emission, the rest are inverse Compton losses by scatter- 
ing of microwave background photons. The amount of en- 
ergy loss of the electrons Qioss ~ 8 • 10'*^ h^Q ergs~^ and 
the power of the flow Qaow onto the relic gives the nec- 
essary minimal efficiency of shock acceleration. The effi- 
ciency needs to be higher than 1% in order to account for 
the radiative energy requirements and for escaping elec- 
trons. This number might be compared to the efficiency 
of shock acceleration in other astrophysical objects: It is 
believed that supernova blast waves have efficiencies of 
1% — 10% for the acceleration of protons as is necessary 
in order to explain the galactic cosmic rays below 10^^ eV 
by supernovae (e.g. Drury et al. 1989; Volk 1997). This 
shows that the assumed cluster efficiency is reasonable. 

2.4. Electron Spectrum 

We apply the theory of plane-parallel shock acceleration, 

because of the large radius of the shock sphere. We mark 
upstream quantities with index 1, and downstream quan- 
tities with index 2. We use a momentum independent 
diffusion coefficient, because of the success of this sim- 
plification in other circumstances (Biermann 1993, 1996; 
Wiebel-Sooth et al. 1997). The distribution function of 
the accelerated electrons cooling by inverse Compton and 
synchrotron emission is 

f{x,p)^Cp"' cxp{-p/p*{x)) (7) 

(Webb et al. 1984), where the normalization C can be ap- 
proximated by being independent of the position x, and 
q = 3R/ {R — 1) is the spectral index of the injection elec- 
tron distribution in three-dimensional momentum space. 
The cutoff momentum is 

X ^ £Po 

^ ^""^ q/ai + {q-3)/a2 + \x\U/{Ka) 
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= : iF+\x\Gy 



(8) 



with U denoting the flow velocity in the inertial frame of 
the shock {Ui — Vs and U2 — Ui/R), k is the diffusion 
coefficient, po is the injection momentum into the accel- 
eration mechanism, and a is the ratio between loss and 
acceleration time scale: 



(9) 



The energy loss time scale rioss mcc/p of an electron with 
momentum p is dominated by synchrotron and inverse- 
Compton losses with: 



4crT / B 



loss 



3m„c V Stt 



(10) 



where Eph is the energy density of the photon field, which is 
dominated by the microwave background. In Sect. 3^ it is 
argued that the main fraction of the synchrotron emitting 
volume belongs to the post-shock region. Thus the electron 
population, integrated over the synchrotron emitting relic 
volume Vol = SD (where S is the surface area, and D is 
the thickness), can be assumed to be 



D 



fip) = S dxf{x,p) 
Jo 
SC 

= TTP'^"^'^ expi-pF) [l-eM'PG2D)] 



(11) 



This spectrum has a break at Pbreak = 1/(^2!?) between 
the spectral index q and and a cutoff at Pcut = ^/F- 

The break in the momentum spectrum at 



T^brcak 



rUc C U2 Tlo 



(12) 



leads to a break in the radio spectrum at t'bioak — 
3eB2 p^,.^^i^/ {2Trmlc^) . Using this relation and solving Eq. 

for the ratio D/U2, which is the time the magnetized 
plasma needed to pass the shock, and therefore the age of 
the relic as a tracer of the shock structure, we get 



tagc — jj — 



D 
T2 



3eB2,cqmcC 



cak 



4 BL ' 



.(13) 



We use a thickness of Z? = 100 ^150 kpc, roughly ^ of 
the width of the projected structure of 1253-1-275, due 
to compression. The age of 1253-1-275 being a shock- 
tracer estimated from the kinematics iagckin = D R/Vg — 
1.5-lO^yr/i^Q'^, and that following from Eq. ^ iage, break > 
4 ■ 10^ yr, using the lowest observation frequency of 151 
MHz as an upper limit to the break frequency inbreak, do 
not fit exactly. Regarding the approximations of this the- 
ory, and the uncertainties of the observational quantities 
only an order of magnitude accordance can be expected. 
But since these two estimates of the relic age depend on 



different observational quantities (iago,kin(fc2obs, ce, D) and 
iagc,brcak(i'broak, ^2,cq)), this is a further test of the theory, 
which could have failed by orders of magnitude. Reasons 
for a discrepancy might also be hidden in neglected prop- 
erties of the acceleration mechanism. For instance if the 
relic has passed the shock already, the diffusion coefficient 
directly behind the shock K2, shock differs from that within 
the relic K2, relic, and the ratio K2,rciic/'«2. shock would enter 
Eqs. |l| and II 

The ratio between break and cutoff momentum of 
12534-275 



c 



Pbrcak 
Pcut 



F 



K.2 



G2D U2D 



q—+q-3 
ai 



(14) 



must be smaller than 0.18, since the radio spectrum does 
not exhibit any break or cutoff over a frequency range 
of about fi orders of magnitude. The term in brackets is 
in the range 1 — 3 and can be neglected for an order of 
magnitude estimate of the diffusion coefficient behind the 
shock: 



K2 



VsDC/R<3- 10^" hso cm^ s"^ 



(15) 



3. Field Compression and Radio Polarization 

3.1. General Considerations 

The high degree of polarization of relic sources should re- 
sult from the compression, which aligns unordered mag- 
netic fields with the shock plane. If the shock processed 
relic is seen at some angle 6 > 0° between the line-of-sight 
and the normal of the shock front, the field structure pro- 
jected onto the plane of sky shows a preferential direction. 
Thus, the resulting radio polarization depends on S and 
R. The direction of the magnetic field should appear to 
be perpendicular to the line connecting the relic and the 
cluster center, because of this compression and projec- 
tion effects. This is the case for 1273-1-275 (Andernach et 
al. 1984). 

In the following we calculate the magnetic field com- 
pression and alignment, and the resulting radio polariza- 
tion of an unordered magnetized plasma blob in an un- 
magnetized gas fiow going through a shock. The (com- 
plex) polarization of synchrotron radiation of an isotropic 
distribution of electrons within a magnetic field seen from 
the z-direction is given as (Burn 1966) 



P 



7+1 By - B1 — 2iBxBy 



7 + 



52 

V 



Bl 



(16) 



The spectral index of the electrons is 7 = 2q: + 1, where 
a is the spectral index of the radio emission. In order 
to average over regions with different field orientations 
and strength, we have to take into account the different 
emissivities 



-sync 



^sync ' 



{Bl+B 



2^- 



(17) 
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Burn has argued that the case a = 1 is a valid simpli- 
fication, which is justified especially in our case, where 
a « 1. Averaging the polarization (denoted by < ... >) 
over a (point symmetric) distribution of fields, weighted 
with the relative emissivity Esync/ < ^sync >, gives the 
observed integral polarization: 



80 



<P> = 



7 + 1 <B[ 



'^>- <Bl 



> 



7+i <B'^> 



<Bl> 



A magnetic flux tube with field strength Bi oriented at g 
some angle 6i to the normal of the shock front will be bent S. 
and amplified by a compression R. From flux conservation 
it follows that 



tan^^2 = R tan 6*1 

B2 = Bi \/cos^0i+R^ sin^ Oi 



Observing this field, which is given by its spherical coor- 
dinates (i?2, ^2; </'2 = 0i)j the viewing angle d results in 



By 



'cos S cos 02 sin 02 — sin S cos O2 

sin <p2 sin 02 | B2 

, sin S cos (/>2 sin 62 + cos S cos 02 



(21) 



Eq. 18 can now be integrated, assuming an isotropic distri- 
bution of the unshocked fields Bi and choosing the com- 
pression ratio of the fields R, which depends on the shock 
compression ratio R and the field orientation. We discuss 
two complementary extreme cases, which give similar po- 
larizations. 

3.2. Weak Fields 

If the magnetic pressure of the relic is small compared to 
the internal gas pressure the compression of the magne- 
tized regions is equal to the compression of the accretion 
shock: R = R. The calculation of the integrals of Eq. 18 
is straightforward if the ensemble average is done over the 
isotropic, unshocked fields, transformed via Eqs. 19, 20, 
and |l| to the shocked fields in the observer's coordinate 
system: 



7 + 1 



7- 



sin^ 6 

fl'i^i ~ sm 



(22) 



Inserting numbers of the relic 1253+275 {5 = 37°, R = 
3.2 ± 0.2, and 7 = 3.36 ± 0.122, derived from a = 
1.18 ± 0.06) gives a polarization of < Pwcak 15%. The 
observed polarization of Pobs — 25% — 30% (Giovannini 
et al. 1991) is higher and might indicate an intrinsically 
ordered field structure of the relic as could be left be- 
hind from ordered field structures within an unshocked 
radio lobe or tail as a progenitor. Or it is due to a larger 
viewing angle oi 6 — 48° — 52°, corresponding to an ac- 
cretion shock radius of rg = 4.0 — 3.6 h'^Q Mpc instead of 




40 

viewing Angle 6 

Fig. 3. Polarization of synchrotron emission of shocked mag- 
netic fields as a function of the viewing angle 5. The strong 
field case is given by solid lines, the weak field case by dashed 
lines. The upper pair of curves (solid & dashed) corresponds to 
a radio spectral index of a = 1, and therefore 7 = 3 and R — 4. 
The lower lines correspond to a = 1.5, 7 = 4 and R = 2.5. Ob- 
served polarizations of radio relics are plotted above the view- 
ing angle predicted from the accretion shock theory (1253-1-257, 
1712+64), simulations of a cluster merger (1706-1-78), or - if 
no data is available - at their best-fit position (0917+75). The 
large projected radius of 0917+75 indicates a large viewing an- 
gle, consistent with the best-fit position. Uncertainties in the 
angles are large. 



Tg = 4.8 /igQ^ Mpc, as given by Eq. 1. Regarding all the 
simplifying assumptions made to derive both estimates 
(Eqs. 1 and ^) of the shock radius, especially the as- 
sumed spherical symmetry, we feel that the derived values 
are quite consistent. 



3.3. Strong Fields 

If the relic is supported by magnetic pressure only, then 
the compression of flux tubes differs from the compression 
of the unmagnetized surrounding medium: The progenitor 
of the cluster relic, an extended radio lobe, was in pres- 
sure equilibrium with its surrounding gas in the upstream 
region. After the radio plasma has passed the shock, it 
expands or contracts rapidly, depending on the field ori- 
entation, in order to achieve pressure equilibrium with its 
environment again, since the downstream flow (in the rest 
frame of the shock) is subsonic. This relates the upstream 
and downstream fields by B^/Bf = P2/P1. Eqs. ^ and 20 
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then give the field compression factor as a function of the 
field orientation: 



I 4R-1 
A-R 



COS^ 6*1 



sm 



The polarization in the strong field case is then 



strong 



>= 



7 + 1 



13-R-7 
R-l 



(23) 



(24) 



Inserting the numbers of 1253+275 gives < Pstrong >= 
16%. This is similar to the weak field case, also demon- 
strated by Fig. ||. Magnetic fields of 1.9 K^q^'^ fJ.G are nec- 
essary within 1253+275 in order to support this struc- 
ture against the thermal pressure at a cluster radius of 
4.8 h'^Q Mpc , assuming a density and temperature as given 
in Sect. 2^ and a magnetic pressure of i3^/(87r). This 
rough estimate depends on an extrapolation of the cen- 
tral density profile of the cluster, and has therefore a 
large error. The magnetic field strength could be lower 
if nonmagnetic pressure such as thermal gas and/or rel- 
ativistic particles also support the radio structure. This 
number is higher than the magnetic field strength de- 
rived from minimum energy arguments of the synchrotron 

2/7 

emitting plasma of 0.5 /I5Q /LtG measured by Giovannini 
et al. (1991). Since these authors use a source depth of 
400ft.^Q^kpc along the hne-of-sight, while we expect only 
125 h'^Q kpc due to the compression of the source by 
the shock, we can correct this number for our model to 
0.7 /igQ^ /LtG. Using the rise of the thermal and (here in the 
strong field case) of the magnetic pressure given by Eq. 
^ we find that the field strength of the unshocked relic 
plasma needs to be 0.18 /xG, if the equipartition value 
gives the relic field strength, or 0.5 /iG if only magnetic 
pressure has to support the relic. Field strengths of the 
pre- and after-shock region should be within these ranges. 

3.4- Depolarization 

A somewhat stronger difference between the polarization 
in the strong and weak field case might arise from inter- 
nal Faraday depolarization, which should only appear in 
the weak field case, since the relic has to be filled with a 
sufficient amount of thermal gas. Internal depolarization 
should be strongest for large viewing angle S, due to the 
larger number of field lines aligned with the line-of-sight. 

We believe that the strong field case is a better de- 
scription of the remnant of a radio galaxy's lobe, since we 
find rough equivalence between the measured equiparti- 
tion field strength and predicted magnetic fields from pres- 
sure equilibrium with the surrounding gas. Therefore, in- 
ternal depolarization should be weak. We do not attempt 
any estimate because of the large uncertainties entering 
the calculation such as coherence scale of the fields, the 
true field strength, and gas density in the relic. 



External Faraday depolarization can lower the ob- 
served polarization if the relic is seen through the central 
magnetized regions of the cluster. The degree of depolar- 
ization would depend on the field strength and reversal 
scale within the intra-cluster medium, where both quan- 
tities have large uncertainties (Kim et al. 1986; Feretti et 
al. 1995; Enfilin et al. 1997; Enfihn & Biermann 1997). In 
order for external depolarization to occur the relic has to 
be on the back side of the cluster, and its viewing angle 6 
must be low. Since then the polarization of the relic would 
be very small anyway, we do not need to correct for this. 

We note that, since the observed polarization of 
1275+273 is high, nearly all of the synchrotron emission 
has to come from the ordered fields of the post-shock re- 
gion as assumed in Sect. 2.4, otherwise the polarization 
would be lower. 

4. X-Ray Emission from Inverse Compton Scatter- 
ing 

The same electrons that emit the synchrotron radiation 
of the relic should inverse-Compton scatter microwave 
background photons to the X-ray range. Goldschmidt & 
Rephaeh (1994) calculated the flux from 1253+275 to be: 



7.1 • 10 



cm^ s ke V V ke V / 



-2.18 



(25) 



Since only a very low thermal X-ray background is 
produced at such a large cluster radius, this inverse- 
Compton flux should be best detectable at keV en- 
ergies. The predicted flux in the ROSAT 0.5-2.4 keV 
band is Fx,o.5-2.4kcV = 1-2 • 10"^ {B2/ ^G)-"^-^^ cuT^ s-\ 
The observed X-ray energy flux from the Coma clus- 
ter is plotted in Fig. 2 of White et al. (1993), and is 
9 • 10^^^ ergcm^^ s~^ deg~^ at the position of the relic 
(assuming a emission temperature of 8 keV). This can 
be translated into an observed flux of i^x,o.5-2.4kcV = 
4 • lO^'^cm^^s^^ coming from the field covered by the 
relic. Comparing this number with the predicted flux im- 
plies a fleld strength of B2 > 0.6 ^G, since most of the 
observed X-ray emission should belong to the X-ray blob 
centered on NGC 4839. This field strength is in good 
agreement with the field strength derived by minimum 
energy arguments and pressure equilibrium, given in Sect. 
3.3. It has to be noted though that if a break in the syn- 
chrotron spectrum below the observed range (150 MHz 
- 4.75 GHz) indicates a break in the relativistic elec- 
tron population, the resulting inverse-Compton fiux in 
the ROSAT band would be largely reduced, since radio 
frequency v translates into inverse Compton energy via 
£^ — 1 keV(i32//^G)~^ (:^/10MHz) in the monochromatic 
approximation of the synchrotron and inverse Compton 
formulae (Blumenthal & Gould 1970). Assuming a break 
in the radio spectral index of 0.5 at 10 MHz or 100 MHz 
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still gives lower limits of B2 > 0.4 /iG or > 0.2 /iG, respec- 
tively. The latter value can be regarded as a hard lower 
limit to the relic field strength. 

1253-1-275 should be an interesting object for future 
X-ray telescopes, such as XMM and AXAF, because its 
large-scale inverse Compton flux might be detected, or 
stronger limits to the magnetic fields and/or the low- 
energy electron spectrum can be derived. The radio halo 
of the Coma cluster should have a larger inverse-Compton 
flux, but there the keV range is polluted by thermal emis- 
sion and has to be avoided by using other energy bands 
(RephaeH et al. 1994, Enfihn & Biermann 1997). 

5. The Origin of Relics 

5.1. Radio Galaxies 

Relics are believed to be remnants of former radio galaxies. 
The problem that the energy loss time scale is frequently 
shorter than the time elapsed since the last resupply with 
energetic electrons from a progenitor radio galaxy is solved 
in our approach due to the acceleration by the large-scale 
accretion shocks. Therefore the progenitor radio galaxy 
needs not to be within the distance from the relic a galaxy 
can travel with a typical cluster velocity within a cooling 
time of the electrons. In order to explain the polariza- 
tion by field compression, it is necessary though that the 
radio plasma from which the present day radio emission 
results was injected within the upstream region. The large 
spatial extent of relics might result from the much lower 
pressure in the upstream regions, so that the radio lobes 
could expand over a larger volume than a similar radio 
galaxy would have occupied in the intra-cluster medium 
inside the accretion shock. 

It is interesting to remember that the minimum mag- 
netic pressure of 1253-1-275 is near equipartition with 
the expected thermal pressure at the shock radius (es- 
timated using an extrapolation of the central density pro- 
file), whereas minimum pressures of radio lobes of active 
radio galaxies in the interior of galaxy clusters usually are 
much lower than the gas pressure at their radii, even if pro- 
jection effects are taken into account (Feretti et al. 1992). 
Several effects can explain this latter pressure imbalance, 
despite the possibility that minimum energies do not de- 
scribe the physical conditions: 

— a filling factor of the observed radio plasma of a few 
percent 

— the ubiquitous presence of relativistic protons within 
the lobes, dominating the energy density of relativistic 
particles 

— the existence of a significant number of low-energy 
electrons radiating below the adopted low-frequency 
cutoff of the minimum energy analysis 

— a significant amount of thermal gas within the lobes. 

Under the hypothesis that relics are old lobes, the filling 
factor and also the low-energy electron population should 



be similar in both types of sources, and therefore not re- 
sponsible for differences of cluster relics and lobes in the 
central regions of clusters. Thermal gas remains to be the 
explanation, and also relativistic protons, because the pos- 
sible rich initial proton population of the lobes might be 
escaped due to the larger age and the higher cosmic ray 
diffusion coefficient of relics compared to lobes inside clus- 
ters. 

Giovannini et al. (1985) propose the Coma cluster 
galaxy IC 3900 to be the origin of 1253-1-275, which should 
have moved to another position on the plane of the sky. We 
point out that the narrow-angle tail (NAT) radio galaxy 
NGC 4789 is a very attractive possible origin, too. First, 
a bridge of radio emission connects it with the central 
ridge of the rehc (Giovannini et al. 1991), indicating a geo- 
metrical connection. Second, 1253-1-275 lies on the (pro- 
jected) line between NGC 4789 and the center of the Coma 
cluster. NGC 4789 is therefore an ideal candidate for be- 
ing the source of magnetized plasma injected upstream 
into the accretion flow. The NAT structure of NGC 4789 
might result from the ram-pressure of the infalling mat- 
ter on the radio jets, enhanced by an ascending move- 
ment of the galaxy. Relativistic electrons from the jets 
should be convected by the flow within the radio bridge 
to the accretion shock, where they are reaccelerated and 
radiate. The line-of-sight velocity of NGC 4789 is higher 
by 1179 km s~^ than Coma's Hubble velocity (Venturi et 
al. 1988). This imphes, if NGC 4789 is really located up- 
stream above 1253-f275, that we see a relic on the back 
side of the accretion shock sphere. This is sketched in 
Fig. ||. If NGC 4789 is on a radial ascending orbit, its 
velocity should be 1480 km/s (assuming 5 ^ 37°), similar 
to that of the infalling matter measured in the clusters rest 
frame Vs = 1400 km/s. The radio plasma, which is reach- 
ing the shock today, was thus injected into the stream 
half a Gyr ago, given that its movement is only due to 
convection by the accretion flow. 

5.2. Large-Scale Fields 

A very different nature of the origin of the magnetic 
fields in the relic is imaginable: from simulations of mag- 
netic field generation and amplification by turbulence in 
the flows of large-scale structure formation field strengths 
of the 0.1 /iC level within cosmological filaments are pre- 
dicted by Biermann et al. (1997). The predicted field 
strengths at cluster accretion shocks are B2 « /xG, being 
in energetic equipartition with the thermal energy den- 
sity. In the inner regions of clusters, field strengths at the 
10 /iG level are expected, as Enfilin et al. (1997) argue to 
be injected from radio galaxies, presuming that their jets 
also contain energetic protons. If field values of B2 ~ fiG 
are typical for the accretion shock region, it has to be ex- 
plained what is special at the positions of relic sources. In 
the case of 1273-1-275 this should be the presence of NGC 
4789 providing preaccelerated electrons injected into the 
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accretion flow, efficiently reacceleratcd at the shock front 
below the galaxy. 

It is remarkable that the 1253+275 complex is in the 
elongation of the Coma cluster, which itself appears to 
align with Coma/Abell 1367 filament of galaxies as de- 
scribed by Fontanelli (1984). It is possible that the po- 
sition of the relic traces the working surface of a large- 
scale stream flowing out of this filament and hitting the 
sea of cluster gas. The observed radio bridge connecting 
1253+275 and the radio halo Coma C (Kim et al. 1989) 
then might be a deposition of magnetic fields from this 
stream, illuminated by relativistic electrons accelerated in 
the shock. 

6. Other Relics 

6.1. General Remarks 

A handful of other cluster relics are known. Their prop- 
erties and that of their host clusters are collected in the 
upper part of Tab. |l|. According to the theory discussed 
in this paper other quantities are calculated and given in 
the lower part of Tab, [l]. Average values for 1253+275 
are also given in Tab. |^but the very recently discovered 
relic 2010-57 (Rottgering et al 1997) is only discussed in 
the text. 

For four relics, namely 0917+75, 1253+27, 1712+64, 
and 1706+78, polarization measurements are available. In 
three cases these polarizations fit well into the predic- 
tion of the accretion shock theory (see Fig. In the 
case of 1706+78 in A2256 the observed polarization is 
much higher than expected. This can be understood if 
the shock is due to an on-going merger event (Fabian & 
Daines 1991), and therefore located at a smaller cluster 
radius. The accretion shock theory predicts a polarization 
of 28% for the relic 2006-56, 10% for 2010-57, 4% for 
1140+203 and 1401-33, and 1% for 0038-096. 

The predicted shock velocity 14, predicted is in all ac- 
cretion shock examples lower than that of the accretion 
shock theory of Kang et al. (1997) by 16% on average. 
This could be due to an assumed temperature drop of a 
factor of less than two from the cluster center to the accre- 
tion shock radius, or due to the need of some recalibration 
of the accretion shock parameters. 

Looking at the distribution of the predicted viewing 
angle 5 one recognizes that large angles are underrepre- 
sented. In an isotropic distribution half of the relics should 
be seen at a viewing angle above 60°. This is clearly a se- 
lection effect, since the sensitive radio surveys of clusters 
with good angular resolution which are needed for the de- 
tection of these sources, usually do not reach projected 
radii of 5 Mpc, as would be necessary in order to include 
the whole accretion shock. We expect therefore that undis- 
covered relics, with strong radio polarizations, are waiting 
in these outer regions of clusters for their exploration. 



6.2. Abell 85: 0038-096 

The spectral index of 0038-096 is not settled in the lit- 
erature: Slee & Siegman (1983) give a — 1.60, Reynolds 
(1986) gives a = 2.90, Slee et al. (1994) 2.37, Joshi et 
al. (1986) 3.03 above a break at 300 MHz, whereas Fer- 
etti & Giovannini (1996) state that a > 1.5. We base our 
estimates mainly on the numbers of the latter authors. It 
is possible that the observed strong steepening is in fact 
a smeared-out cutoff, resulting from the limited accelera- 
tion power of a weak accretion shock. Since the projected 
relic position is close to the cluster center the viewing an- 
gle is close to zero and therefore no visible polarization 
resulting from field compression is expected. The nec- 
essary radio efficiency Tjradio of the shock is > 1% and 
thus reasonable. The predicted magnetic field strength 
-B2, predicted = 2.6 /xG /igQ^^^ from assumiug pressure equi- 
librium between the post-shock gas and the relic is again 
higher than the equipartition (or minimum energy) value 

2/7 

given by Feretti & Giovannini (1996) -B2,oq < 1.4/i G/igQ 
(corrected by i?^/^). Lima Ncto et al. (1997) report the 
appearance of an X-ray blob in the X-ray image of A85 
centered on 0038-096, possibly inverse- Compton scattered 
microwave background photons. If so, the magnetic field 
strength could be derived directly from an estimate of this 
X-ray excess. 

6.3. Abell 786: 0917+75 

The megaparsec cluster radio relic 0917+75 is located 5 
Mpc h'^Q away from the center of A786 within the plane of 
sky. Unfortunately no temperature or velocity dispersion 
of A786 is available in the literature. Thus, the expected 
shock radius and therefore the polarization cannot be pre- 
dicted. However the very peripheral position indicates a 
large viewing angle, consistent with the high degree of 
the observed polarization of 48 — 60%. A viewing angle of 
72° — 63° corresponding to an accretion shock radius of 
5.6 — 5.3Mpc/i^Q^ would be in agreement with this polar- 
ization. Such a large shock radius indicates a deep grav- 
itational potential. The corresponding virial temperature 
13 — 15 keV (Eq. 1) is in the upper range of observed tem- 
peratures of clusters, and therefore possible. However the 
shock structure might in fact be also due to an accretion 
flow, which is centered on and stops at the boundary of 
the supercluster Rood #27, which A786 is a member of. 
This could explain the large projected shock radius. 

The direction of the projected magnetic field observed 
by radio polarization is mainly perpendicular to the axis 
connecting the center of A786 and 0917+75, as it should 
be if the field is compressed at a spherical shock centered 
on the cluster. But a weakly visible spiral structure in the 
projected fields could result from an intrinsic order which 
the magnetic fields had before they passed the shock. 

The relic 0917+75 exhibits a break in its spectrum 
(Harris et al. 1993). Below 300 MHz the spectral index 
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Table 1. The upper part of the table is a collection of measurements of properties of cluster radio relics and their host clusters. 
The references are indicated in the third column in the sequence of the eight examples. The lower part gives estimates, using 
the above properties and the formulae of this paper, as indicated in the third column. A discussion of the individual sources 
and various assumptions made in order to get the given values can be found in Sect. ^ 



Observed Quantities: z: redshift of the cluster. kTohs'- observed central gas temperature. no,o, ''core, P'- parameters of the 
/3-model, describing the radial electron density profile of the cluster no = ne,o [1 + {r /rcoic)^]~^'^^^ ■ '"projected: radius of the 
position of relic, projected into the plane of sky. a: spectral index of relic's radio emission. P: average observed polarization. 
B2,eq'- equipartition field strength of the relic, the index 2 refers to the post-shock region, where to the visible part of the 
relic belongs. The values of the literature are changed due to the effect of a smaller volume than assumed there according to 
S2,eq ~ R^^'^ , where R is the compression ratio. D: thickness, roughly estimated by dividing the observed projected thickness 
by the compression R. S: surface area of the relic. Since in all except one cases the relic structure can not be deprojected, we 
use the observed projected surface as a rough estimate, f break: break frequency of relic's radio spectrum. (: ratio of break- to 
cutoff momentum of the electron population, estimated from the ratio of break- and cutoff frequency of the radio spectrum by 
taking the square root. Qradio: relic's radio power, usually assuming a single power law between 10 MHz and 10 GHz. 
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Estimated Quantities: Vs'. shock radius. Vs'. shock velocity. 5: viewing angle given by the estimate of the shock radius Va and 
the projected radius of the relic rproj. R'- shock compression ratio, estimated from the spectral index. P2/P1, T2/T-1: pressure and 
temperature jump across the shock. nc,i: upstream side electron density, derived from the /3-profile and the compression ratio. 
kT\: temperature of the infalling gas, derived from the post-shock temperature kT2, which is assumed to be half of the central 



temperature fcTobs (see Sect. 2.2). predicted: shock velocity, predicted from the post-shock temperature and the compression 
ratio. -Bi, predicted, predicted : Up- and dowustream field strength which would be in pressure equilibrium with the surrounding 
gas. iJi.eq: strength of the magnetic field of the progenitor of the relic source, if the magnetic pressure was increased by P2/P1 
to -B2,eq in the shock. Qnov- kinetic power of the flow onto the surface of the relic. ?7radio: radio efficiency Qradio/Qflow tagc,i<in: 
age of the relic being behind the shock from D/U2-, where U2 = Vs/R is the post-shock velocity, tago, break: age of the relic 
being behind the shock from D/U2, where this ratio is derived from the radio break frequency. K2: diffusion coefficient of the 
post-shock region. Pstrong: polarization expected in the strong field case, for a viewing angle 5. (5strong, J's, strong: viewing angle 
and the shock radius computed to be consistent with the observed polarization P. 

Remarks: 

The lower limit to the spectral index translates into limits of dependent quantities. But the measured steep spectral index 
might also be due to the appearance of the high frequency cutoff, so that these limits become estimates. ^'^^ A spectral index 
smaller than unity indicates, that the range observed in the radio is below the break frequency. The spectral index above the 
break is steeper by 0.5. A after-break spectral index of a = 1.3 was used in these formulae. No lower limit to the ratio (" between 
cut and break frequency and therefore no diffusion coefficient can be estimated in this case. The results of the accretion 
shock theory are — 4.9Mpc /i^Tg^, Vs = 2280 km/s, and Pstrong ~ 0.78%. The latter is excluded by the observed polarization 
of P = 20%. The accretion shock model therefore fails; but as discussed in Sect. the shock should result from a merger. We 
therefore inserted the best-fit numbers from a hydrodynamic simulation of such an event, is set to rs, strong, and 5 axe taken 
from Roettiger et al. (1995). Since this shock is in the interior of the cluster, kT2 = fcTobs was used. For this relic a 

deprojection of the geometry was possible as described in the text, and we therefore use these values. '■^^ The velocity dispersion 
given in Fadda et al. (1996) is translated into a temperature with the help of the dispersion-temperature relation given in Bird 
et al. (1995). 

References: |: Briel et al. (1991), §: Briel et al. (1992), |: Burns et al. (1995), |d[ David et al. (1993), |[ David et al. (1995), 
I Fadda et al. (1996), |: Feretti & Giovannini (1996), R_Feretti et al. (1997), |_Gavazzi (1978), |: Gavazzi & Trinchieri (1983 
|: Gavazzi & Jaffe (1987), | Giovannini et al. (1991), m[ Goss et al._(1982),^ Goss et al. (1987), § Harris et al. (1993), [ 
Jones & Forman (1984), H: Joshi et al. (1986), ^ Knopp et al. (1996), |: Lima Neto et al. (1997), |: Rottgering et al. (1994), [ 
Rottgering et al. (1997), R Zabludoff (1993) 



is 0.6 and above it is steeper by the canonical value 0.5, 
as is expected from a steepening of the spectral index of 



the electron population by 1, predicted in Sect. 2.4. The 
corresponding compression ratio of i? = 3.5 is close to the 
maximum possible value i? = 4 of a strong nonrelativistic 
shock. 

The projected peripheral position of 0917-1-75 al- 
lows a sensitive search for inverse-Compton scattered mi- 
crowave background photons in the X-ray range. Harris 
et al. (1995) observed the relic with the PSPC instrument 
on the ROSAT satellite and got an upper limit to the 
flux in the 0.5-2.0 keV band of 2.3 • lO"!'' ergcm^^ g-i 

This implies that the magnetic field strength must be 
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B2 > 0.8 ^G, shghtly larger than i?2,cq = O-^^Gh^Q de- 
rived from minimum energy considerations and various 
assumptio ns, which were discussed in Harris et al. (1995) 
and Sect. 



5.1 



of this article. 



6.4. Ahell 1367: 1140+203 

The diffuse radio source found in A1367 is usually called a 
halo, but it is more probably a cluster relic because of its 
nonccntral location and asymmetric shape. The distance 
to the center of the galaxy distribution is 0.8 Mpc /i^q^ 
(Gavazzi 1978), but to that of the X-ray emission it is 
l.lMpc/i^Q^ (see map in Gavazzi et al. 1995). We use the 



latter distance, since the X-ray emission should be a better 
tracer of the gravitational potential. The X-ray emission 
is elongated into the direction of the relic, possibly trac- 
ing the influence of the large-scale structure potential, or 
the main direction of the accretion flow onto the cluster, 
or both. The radio relic could be contaminated or related 
to three close irregular galaxies, which have radio trails 
behind them (Gavazzi & Jaffe 1987), indicating that they 
are falling inwards. Strong on-going star formation, a re- 
cent supernova in one of them, and young and abundant 
H H regions, which are aligned as if they were formed 
by bow shocks (Gavazzi et al. 1995), could be triggered 
by the sudden change of the environment of these galax- 
ies during the passage of the accretion shock. A starburst 
might increase the accretion within these galaxies onto a 
possible central black hole (similar to the argumentation 
of Wang & Biermann 1998) and therefore can trigger the 
ejection of radio plasma. 

The measurements of the spectral index and radio flux 
of the relic differ in the literature (compare Gavazzi 1987, 
Hanisch 1980, and Gavazzi & Trinchieri 1983) and there- 
fore any derived numbers have to be used with care. Espe- 
cially the radio power Qiadio is an extrapolation of a very 
steep spectrum (a — 1.9) down to 10 MHz and proba- 
bly an overestimate. Due to the low expected viewing an- 
gle only a small polarization of 4% is expected. Gavazzi 
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& Trinchicri (1983) detected some X-ray emission from 
the region of the rehc. Assuming this to result from in- 
verse Compton emission they estimate a magnetic field 
strength of 1.2 fiG, of the same order as the equiparti- 
tion value of -B2,cq = 2 /xG h^Q and the predicted strength 
of i32, predicted = 2.4 //G ft-gp^ for pressure equilibrium with 
the surrounding. 

6.5. Abell 2255: 1712+64 

The cluster A2255 has a number of similarities to the 
Coma cluster: there is evidence that the cluster is in a state 
of merging (Burns et al. 1995), a radio halo is present, and 
also a peripheral radio relic without any obvious parent 
galaxy (Feretti et al. 1997). The spectral index of the relic 
is fairly constant along the structure. The upper limits 
to the polarization of 9% at 20 cm and 2% at 90 cm do 
not contradict the predicted polarization of 3% seriously, 
taking into account e.g. the rough assumptions which en- 
tered the estimate of the shock radius, and the possibility 
that external Faraday depolarization has lowered the de- 
gree of the radio polarization if the relic is seen through 
the denser magnetized central intra-cluster medium. This 
relic therefore fits also well into the accretion shock theory. 

6.6. Abell 2256: 1706+ 78 

The radio emission reveals a very complex situation within 
A2256: a mini- halo is visible in the map of Bridle & Fo- 
malont (1976), and several head-tail radio galaxies can 
be found. One NAT source (in the notation of Bridle & 
Fomalont (1976): source C) is one of the longest head- 
tail sources known. The tail extends over 700 kpc /i^g^ , is 
slightly bent and exhibits a kink 570 kpc Ii^q behind the 
head (Rottgering 1994). The bending can be understood 
as the gravitational influence of the cluster mass to the 
trajectory of a fast radio moving radio galaxy, at which 
the tail traces the path. The kink might result from a sud- 
den change of the surrounding density or velocity field of 
the background gas (Rottgering et al. 1994), as it might 
happen if the radio galaxy passes a shock. The projected 
position of the kink is within a very extended irregular 
and sharp-edged region of diffuse radio emission in the 
north-west region of the cluster. This region consists of 
two extended sources (G and H in the notation of Bri- 
dle & Fomalont (1976)), which might be two independent 
relic sources with similar properties. 

The spectral index of 0.8 of these possible cluster rehcs 
indicates that the injection electron spectrum below the 
break in the momentum distribution are observed. The 
spectral index of the synchrotron emission above the break 
should be a — 1.3. The break momentum (Eq. |l^) is 
highest for a large shock velocity Vg and a thin magne- 
tized post-shock region (small D). The latter could ex- 
plain why the emission region is relatively sharp-edged. 
If the thickness D is much smaller than we assumed, the 



kinematic age of the relic iage,kin would be much smaller 
than 10 yr, in accordance with iage, break 

< 1.4 • 10** yr. 

The polarization of 20% is much too high for the expected 
viewing angle S — 9° oi a relic located at the canonical 
radius of the accretion shock — 5 Mpc h'^Q , which pre- 
dicts only 1% polarization. A consistent viewing angle of 
^strong ~ 43° would imply that if the normal of the shock 
front is pointing to the cluster center, the shock radius is 
Ts = 1 Mpc h'^Q . This would be in agreement with the aver- 
age direction of the polarization, which is oriented roughly 
perpendicular to this assumed normal of the shock front. 
A shock in the interior of the cluster can only be explained 
by a merger event. Independent evidence for such a colli- 
sion can be found in the X-ray morphology, in the com- 
plex temperature structure and in the elongated galaxy 
distribution. An extensive discussion of the observational 
signs reported in the literature is given by Roettiger et 
al. (1995). They compare these observations with simula- 
tions of cluster mergers and conclude that the entire data 
are best reproduced by a merger of two clusters with mass 
ratio 1:2 and a relative velocity of 3000 km/s. The sepa- 
ration of the cluster centers is 0.75 Mpch^Q. The smaller 
cluster is infalling from the north-west and is approaching 
the observer. The angle between the line-of-sight and the 
merger axis is in the best-fit model in the range of 40° to 
45°, in agreement with ^strong = 43° estimated from polar- 
ization, without taking any possible influence of external 
Faraday depolarization into account. 

The temperature fcTobs given in Tab. |l| is an average 
value of a complex distribution, biased by the heat pro- 
duction of the shock. Since the infalling subcluster has a 
much higher temperature {kTi « 4 keV; kT2 = fcTobs is 
assumed here) than the accretion streams hitting all other 
known cluster relics {kTi « 1 keV), the moderate shock 
compression ratio R = 2.9 at this shock, which has the 
highest shock velocity of our sample, can be understood 
to result from the high pressure of the upstream matter. 

The predicted shock velocity derived from the up- 
stream temperature and pressure jump is Ks, predicted = 
3050 km/s, being in agreement, and completely indepen- 
dent of the above mentioned merger velocity from the 
hydrodynamic simulations. The merger velocity corre- 
sponds to the infall velocity measured in the cluster rest 
frame Vs, which is usually lower by {R — 1)/R than the 
shock velocity 14 . But in case of a massive merger it is un- 
realistic to believe that A2256 behaves like a single body, 
the downstream flow does surely not rest in the cluster 
reference frame. We therefore use a realistic velocity of 
V; = V; = 3000 km/s. 

The line-of-sight velocities of the head-tail radio galax- 
ies also point towards the observer, and especially the ob- 
viously fast moving source C seems to originate from the 
infalling cluster. Rottgering et al. (1994) have fltted tra- 
jectories, calculated for a gravitational potential derived 
from X-ray data to the tail of source C, and concluded 
that the initial velocity had to be in the range of 2000-3500 
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Fig. 4. Central cluster temperatures vs. spectral index of the 
radio relics. The dotted lines are temperature-spectral index 
relations, computed from the temperature jump given by Eq. 
^ which for a given temperature of the infalling gas kTi (and 
assuming kT2 — |fcTobs) is a function of the spectral index a 
only. The temperature of the infalling gas fcTi is indicated for 
the curves. 



km/s, in accordance with the value of the relative cluster 
velocity given by Roettiger et al. (1995). From the undis- 
turbed straightness of the source Rottgering et al. (1994) 
derive that the turbulent motion of the cluster medium 
has a velocity < 200 km/s, much too low to explain the 
diffuse emission resulting from large-scale turbulent mo- 
tion. However, this is not in contradiction to the merger 
shock model for this relic, where violent velocity changes 
occur only in the thin shock region, which might have dis- 
torted the tail of source C at the position of the kink. We 
note, that a brightening of the radio tail is visible there, 
as expected to result from the action of the shock onto the 
tail (Rottgering et al. 1994). Source C should have moved 
from the kink within 2 — 3 ■ 10^ years to its present day po- 
sition. This is comparable or higher than a typical cooling 
time of radio emitting electron, and explains the observed 
fading of the tail. 

The details of the complex scenario might differ from 
the above picture in that the shock might be oblique 
(Roettiger et al. 1995), or the impact parameter of the 
merging clusters is non-zero (Briel et al. 1991). 

6.7. Abell 3667: 2006-56 and 2010-57 

The giant cluster radio relic 2006-56 is peripherally lo- 
cated within the direction of the elongation of A3667. This 
elongation is seen in X-rays, but also contains a subgroup 
of galaxies (Knopp et al. 1996). Because of the low X-ray 



background, 2006-56 would be an ideal candidate for a 
future sensitive search for inverse Compton flux (by, e.g., 
AXAF, XMM). Goss et al. (1982) separated two extended 
radio sources at this position, located close to each other: 
2006-56 and a diffuse ridge. A recent sensitive observation 
by Rottgering et al. (1997) shows that there is in fact a sin- 
gle Z-shaped structure, extending over 2.6 Mpc /i^g^, with 
constant spectral index along the main axis of the relic. 
The Z-shape might be the structure of the radio plasma, 
since we do not see the relic edge-on. The side more distant 
from the cluster center is sharply edged and has a flatter 
spectral index, as is expect ed t o happen on the upstream 
side as explained in Sect. ^.l[ . The spectral index drops 
from a = 1.1 in the central region of the relic to « 0.5 
at this edge. This indicates strongly the presence of par- 
ticle acceleration, so that we can use the thickness of the 
rim of lower spectral index in order to estimate the thick- 
ness of the relic. The projected thickness derived from the 
rim is approximately 100 — 200kpc/i^Q^. The region with 
steeper spectral index on the inner side has also a pro- 
jected thickness of « 200 kpc /i^g^. Assuming a viewing an- 
gle oi 5 — 48°, deprojection gives D « 150 — SOOkpc/i^Q^. 
The projected area of the relic is 1.6 Mpc^ ■ ^ stripe of 
150 kpc X 2.6 Mpc belongs to one edge. Deprojecting the 
remaining area gives a relic surface of 1.8 Mpc^ /i^q^. The 
infalling gas seems to be relatively cool, with kTi « 0.4 
keV, but this number has a large error since the tempera- 
ture jump is very sensitive to small variations in _R, since 
R is close to its maximum value (see Eq. ^ or Fig. 

Rottgering et al. (1997) report a pressure imbalance 
of one to two orders of magnitude, between the gas pres- 
sure at the projected cluster radius of the relic and its 
equipartition pressure. This discrepancy should vanish if 
the lower gas pressure at the shock radius and the higher 
radio plasma pressure of a compressed, flattened structure 
is taken into account. Therefore we flnd good agreement 
between equipartition field strength and predicted field 
strength. 

Surprisingly, the radio image of Rottgering et al. 
(1997) reveals a second cluster radio relic on the other 
side of the cluster (2010-57). These authors mention that 
both relics might be the two lobes of a former, now inac- 
tive, radio galaxy, possibly the cD galaxy. This must have 
been gigantic, since the projected distance is 5.2 Mpc h'^^ ■ 
The expected diameter of the accretion shock sphere of 
8 Mpc h'^Q would be a lower limit to the extent of this ra- 
dio galaxy in our scenario, since nowadays we would see 
a backflow of the radio lobe from larger distances. In or- 
der to achive such an extent, the radio galaxy must have 
been very powerful and the ambient medium more tenu- 
ous than today. If both structures are really produced by 
one galaxy, this should have happened during the early 
period of violent quasar activity, when the intra-cluster 
medium was thinner. 

If the two relics are independent, in the sense that 
their radio plasma is from different sources, but their po- 
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sitions are correlated due to the membership of A3667 
to a large-scale filament, then the head-tail radio galaxy 
B2007-569 can be the source of the shocked radio plasma 
in 2006-56. It is located between the cluster center and the 
relic and is moving towards us with a speed of 1200 km/s 
{z = 0.05257; Sodree et al. 1992). If it is on a radially 
falling orbit, which traversed the shock sphere at the lo- 
cation of the present day relic, it should have a velocity of 
1800 km/s (using S = 48°), higher than the expected infall 
velocity of matter ( measured in the cluster rest frame) at 
the higher radius of the shock. Travelling with this veloc- 
ity it should have passed the shock a Gyr ago. This is an 
order of magnitude longer than a typical cooling time of 
radio emitting electrons and therefore the expected trace 
of radio plasma on the path of that galaxy became in- 
visible. It is amazing that Rottgering et al. (1997) found 
marginal evidence for a bridge between this galaxy and 
the relic, which might indicate reacceleration. 

We predict a polarization in 2006-56 of 28% from our 
accretion shock model, and 10% in 2010-57 due to its 
smaller expected viewing angle of 30°, assuming the same 
spectral index and shock compression ratio as for 2006- 
56. The observation of these polarizations would be an 
independent test of the theory, and therefore could verify 
general assumptions made about large-scale flows in the 
Universe. 

6.8. Abell S753: 1401-33 

S753 is a poor cluster. Its low velocity dispersion of 536 
km/s (Fadda et al. 1996) indicates a low central gas tem- 
perature of fcTobs = 2.5 keV, applying the dispersion- 
temperature relation given in Bird et al. (1995), and there- 
fore a flat gravitational potential compared to the other 
clusters. The spectral index a = 1.4 of the relic 1401- 
33 corresponds to a shock compression ratio of R — 2.7. 
This relatively strong shock at this poor cluster can be 
understood if the accreting matter has a temperature of 
kTi « 0.5 keV. This is lower than what seems to be typ- 
ical for the other clusters (fcTi w 1 keV, see Fig. |). The 
lower preheating of the flow onto this cluster might in- 
dicate a poorer cosmological environment. Since the sky 
position of the relic 1401-33 is close to the center of S753, 
the viewing angle should be small and only 4% polariza- 
tion is predicted. 

7. Conclusions 

From simulations of the structure formation of the Uni- 
verse, large accretion shocks around clusters of galaxies are 
expected. If magnetic fields are present at these locations, 
particle acceleration should take place. Protons might be 
accelerated to the highest observed energies, as demon- 
strated by Kang et al. (1997), because of their weak energy 
losses and their possible injection out of a thermal popula- 
tion. If a population of relativistic low-energy electrons is 



also present, they would be accelerated, too. High-energy 
ultra-relativistic electrons lose energy by inverse-Compton 
scattering and synchrotron emission. This allows to see 
them in the radio or X-ray regime and therefore verify the 
existence of the large-scale formation shock structures. 

At peripheral positions with respect to the cluster cen- 
ter regions of diffuse radio emission are found, the so-called 
cluster radio relics ( Figs. |l| and They are assumed to 
be remnants of radio galaxies. Their energy supply was a 
long lasting outstanding problem, since their electron cool- 
ing times are usually too short to allow any of the nearby 
possible former parent galaxies to have moved from the 
relic to its present position within that time. Electron ac- 
celeration has to take place within them, but turbulent 
acceleration from galactic wakes fails to explain this, due 
to the low galaxy densities in these regions. 

If these relics are close to the accretion shock, a shock 
efficiency of only 0.1 — 5% is sufficient to power them. The 
accretion shock itself should collect relics since every vol- 
ume filled with radio plasma, being injected into the inter 
galactic medium above such a shock by a radio galaxy, 
would be dragged into the shock by the fast motion of the 
infalling gas. After passing the shock, the relic remains 
close to it for some time, since the growth rate of the 
shock radius is much slower than the accretion velocity. 
Buoyant motion of the light radio plasma embedded in 
the heavier surrounding gas might help to keep it longer 
there. 

The relic is proposed to be radio illuminated by rela- 
tivistic electrons (re) accelerated at the shock and there- 
fore tracing its position. The radio spectrum of such relics 
is steeper than expected for electron acceleration within 
strong shocks, indicating that the temperature of the in- 
falling matter should be « 0.5 - 1 keV (see Fig. ^. Thus, 
preheating of this gas should have taken place. The ac- 
creting gas has a typical density of ne,i ~ 10~^ cm~^ ^5o^- 
1706-1-78 in A2256 of course does not fit into this temper- 
ature scheme, due to its different physical nature, being 
located at a merger shock front in the interior of the clus- 
ter. Matter which has fallen onto the large-scale cosmic 
sheets should be shock- heated to about 0.1 keV and have 
a density of < 10^^ cm~'^ far away of clusters. Adiabatic 
compression and internal shocks of the converging flow 
within these sheets onto clusters of galaxies should heat 
the gas further to the above temperatures of 0.5 — 1 keV. 

During the passage through the shock, the magnetic 
flelds are amplified and aligned with the shock plane. The 
radio emission of a relic should therefore exhibit polariza- 
tion properties depending on the viewing angle. Assuming 
a spherical accretion shock with a radius depending on the 
depth of the cluster potential, the viewing angle and the 
resulting polarization can be estimated (Fig. ^). These 
agree well with measured polarizations for relics in two 
clusters: 1253+275 in Coma, and 1712-1-64 in A2255. For 
0917-1-75 in A786, the large projected radius would imply 
a deep gravitational potential and therefore a very high 
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cluster temperature, which is not yet measured. 0917+75 
might also trace an accretion shock of the stream onto the 
supercluster Rood #27, which A786 is a member of. This 
could also explain the large apparent radius of the relic. 

In case of 1706+78 in A2256 the measured polarization 
exceeds the expected one significantly, indicating that the 
viewing angle is much higher than if the relic were lo- 
cated at the accretion shock. A consistent shock radius 
would be 1 Mpc h'^^ , which means that an internal shock 
of the cluster powers the relic, which should result from 
a highly developed merger event. Roettiger et al. (1995) 
discuss several apparent signs of such an event in A2256, 
and fit simulations to the X-ray data. Their best-fit value 
of the angle between the merger axis and the line-of-sight 
within our theory of field compression predicts exactly the 
observed polarization. 

The directions of the projected magnetic fields are 
roughly perpendicular in all above examples to the di- 
rection pointing to the cluster center, as is expected. 

For the other five relics of our sample no polariza- 
tion measurements are available. We predict « 30% po- 
larization for 2006-56 in A3667, « 10% for 2010-57 in 
A3667, and less than 5% for 0038-096 in A85, 1140+203 
in A1367, and 1401-33 in S753. Polarization measurements 
of these sources could verify the existence of cluster accre- 
tion shocks and help us to reveal an important piece of 
the cosmological structure formation puzzle. 

The time elapsed since the magnetized plasma has 
started to pass the shock can be calculated in two in- 
dependent ways. Both estimates give ages of the order of 
some 10® yr, and are therefore a further justification of 
our model. 

Since our sample contains mostly relics seen under 
viewing angles below 60°, we expect that more relics might 
be discovered in the outer regions of clusters. Any radio 
search within these large areas should be guided by the 
expected correlation of the large-scale structure filaments 
and typical inflow directions into clusters. 

We conclude that the properties of cluster radio relics 
are naturally explained if they are understood to trace 
some of the giant shock fronts of the cosmological large- 
scale motion of the on-going structure formation. This 
demonstrates for the first time the existence of these theo- 
retically predicted shocks. Estimated magnetic fields, tem- 
peratures and densities of the accreting matter fit into the 
structure formation scenario. 
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